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Abstract 

The nucleosynthesis and production of radioactive elements in SN 1987A are re¬ 
viewed. Different methods for estimating the masses of 56 Ni, 57 Ni, and 4i Ti are 
discussed, and we conclude that broad band photometry in combination with time- 
dependent models for the light curve gives the most reliable estimates. 
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1 Introduction 


Two of the outstanding issues for SNe are the nucleosynthesis and the explo¬ 
sion mechanism. Although SNe have been known from the 1950’s to be the 
most important sources of heavy elements, there is little quantitative evidence 
for this. SN 1987A has in this respect been a unique source of information 
thanks to the possibility to obtain spectral information at very late epochs 
when the ejecta is transparent and the central regions observable. The nature 
of the explosion mechanism is still to a large extent unknown, and constraints 
from observations are badly needed. The main such constraints are provided 
by the hydrodynamic structure, e.g., the extent of mixing, and the masses of 
the iron peak elements, in particular the radioactive isotopes (e.g., [Timmes! 


Woosley, Hartmann, fe Hoffman, 1 lJTJ(3]: [Kumagai et ah, 199cS|; [lxifonidis, Plewa, 


Janka, fc Miillcr, 2000|) . 


2 Abundances 


To model the spectra and light curve a complicated chain of thermalization 
processes from the gamma-rays and positrons in the keV - MeV range, to the 
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Table 1 

Main radioactive decays in SN 1987A 


Decay 



Time scale 

Epoch when dominating 

56 N f 

-> 56 Co + 7 

56 Co 

—► 56 Fe + 7 

-*• 56 Fe + e+ 

8.8 d 

111.3 d 

0 - 18 d 

18 - 1100 d 

57 Ni 

-»• 57 Co + 7 

57 Co 

-*• 57 Fe + 7 

2.17 d 

390 d 

1100 - 1800 d 

44 Xi 

-> 44 Sc+ 7 

44 Sc 

-»• 44 Ca + 7 

87 yrs 

5.4 h 

1800 d 


44 Ca + e+ 


observed optical and IR photons has to be understood. The details of this are 
discussed in a series of papers, ([Kozina fe Fransson, 1992j , |1998a| Jb|). See also 
Fransson (1994| ) and |dc Kool, Li, fc McCray (1998|) . Here we only summarize 
the main points. 


After a couple of days the main energy input to the SN ejecta comes from ra¬ 
dioactive decay, first 56 Ni followed by 56 Co. Later than ~ 1100 days 57 Co takes 
over, while at very late epochs, >2000 days 41 Ti decay becomes dominant. Ta¬ 
ble 1 summarizes the different decay routes, the emitted particle (gamma-rays 
or positrons), the exponential decay time scale and the approximate epoch 
when the various decays dominated the energy input to SN 1987A. 


The gamma-rays produced by the radioactive decay loose their energy by 
Compton scattering off bound electrons, producing a population of non-thermal 
electrons in the keV range. These, as well as the positrons, are slowed down 
by scattering off thermal, free electrons, ionizations and excitations. This cas¬ 
cade has to be calculated either by Monte Carlo or by solving the Boltzmann 
equation, as discussed in |Kozma fe kransson (1992] ). Knowing the heating and 
ionization rates, one can then calculate the temperature and ionization state 
of the gas, and from this the emission from the different nuclear burning zones. 
Finally, by solving the radiative transfer one can calculate a spectrum, which 
can be compared to the observations. 


In figure 1 we show as an example the light curve of the Ha line compared 
to data by Danziger et al.(1991) and the SINS/HST collaboration (Wang et 
al.1996; Chugai et al.1997). For the explosion models we have tested both the 
11E1 and 14E1 models by Shigeyama & Nomoto(1990) and the 10H model by 
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Fig. 1. Left panel: Evolution of Ha from Danziger et al. and from SINS/HST 
observations together with model calculations. The flat tail is mainly a result of 
freeze-out rather than input from 4i Ti, as can be seen from the dominance of the 
low density envelope (dashed line), compared to the higher density core (dotted). 
Right panel: Evolution of the ’6300 A’ line. At early epochs the [O I] AA6300, 6364 
lines dominate (dotted line), while most of the flux at epochs i>1400 days originates 
from an Fe I line at 6300 A. The oxygen mass was in this model 1.9 M 0 . 

Woosley (1988). Thge difference between these models are, however, marginal. 
The agreement between the model and the observations is very good over the 
whole evolution, which gives confidence in the whole procedure. A problem 
here is that a substantial fraction of the hydrogen may be at high velocity, 
where it is difficult to detect in the faint extended line wings. An accurate 
reproduction of the line profile is therefore a must in this case. In our models 
we find that the total mass of hydrogen rich gas is 7.7 ± 2 M@, of which 
~ 2.2 M 0 is mixed within 2000 km s -1 . Of the hydrogen rich gas ~ 3.9 M 0 
is helium, while an equal mass is hydrogen. 

The mass of the helium dominated zone is consistent with ~ 1.9 M 0 , of which 
~ 1.8 Mq is pure helium, with a further ~ 3.9 M 0 of helium residing in the 
hydrogen component. A possible caveat is again that a substantial mass could 
be hidden at velocities > 3000 km s _1 . 

Oxygen is important both as a diagnostic of the evolutionary calculations, the 
progenitor mass, and because of its large abundance. In figure 1, right panel, 
we show the corresponding light curve for this line. A technical point is here 
that while the feature at ~ 6300 Ais dominated by the [O I] AA6300, 6364 
doublet at epochs <1000 days, it is at late time instead dominated by an Fe I 
line at 6300 A. This illustrates the importance of a self-consistent calculation 
of the spectrum, taking all important ions and transitions into account. The 
total oxygen mass we find is ~ 1.4 M 0 , but the permissible range is probably 
as large as 0.7 — 3.0 M 0 . One can from this with good confidence rule out 
progenitor models with ZAMS mass less than 15 M 0 , because of the low 
oxygen mass (<0.5 M 0 ). 

In table 2 we summarize the most important masses derived from the spectra 
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Table 2 

Masses of some important elements and isotopes in SN 1987A. 


Isotope 

Mass ( Mq) 

Method 

Note 

1 H 


3.9 

Spectrum 


4 He 


5.8 

Spectrum 


16 O 


1.9 

Spectrum 


56 Ni 

(56 Fe ) 

0.069 

Bol. l.c. 

±0.003 Mq 

57 Ni 

( 57 Fe) 

0.003 

UBVRI l.c. 


£ 

00 

LO 


0.006 

Spectrum 


44 Ti 

( 44 Ca) 

0.0001 

UBVRI l.c. 

(0.5 - 2) x 10 -4 Mq 


and light curves (bolometric and broad band UBVRI). The latter are discussed 
in section 3. Other elements are discussed in IKozma fe Fransson (1998b|) . 


3 The radioactive isotopes and the light curve. 


To estimate the masses of the three most abundant radioactive isotopes formed 
in the explosion, 56 Ni, 57 Ni and 4 i Ti, one can use several approaches, each with 
its special advantages and disadvantages (Fransson & Kozina, in preparation). 
The best known, and for the 56 Ni mass most reliable, is the bolometric light 
curve. As shown by e.g., |Bouchet et al. (1991|) , this yields a 56 Ni mass of 
0.069 M e . 

For 57 Ni, or rather, 57 Co, which dominates the light curve between 1100 - 1800 
days, this is not so straightforward, and a simple approach similar to that of 
56 Co, yielded a 57 Ni mass at least twice as large as other methods, and much 
higher than could be produced by the nucleosynthesis models flSuntzetf et 
al., 1992| ; |Bouchct, Danzigcr, fc Lucy, 199 1|) . This was shown by [FranssonHc 
Kozina, (1993| ) to be a result of time dependent effects in the light curve. Be¬ 
cause of the decreasing density and ionization, the recombination time scale 
(oc l/(a rec n e )) as well as the cooling time scale, become longer than the ex¬ 
pansion time scale, t. The gas is therefore not able to recombine and cool at 
the same rate as the radioactive input takes place, and only a fraction of the 
radioactive input is released as radiation. The difference is released later as 
recombination radiation, giving a higher luminosity than the instantaneous 
input. The bolometric light curve therefore gives an overestimate of the ra¬ 
dioactive mass if a simple steady state calculation is employed. 


In figure 2 we show a time dependent calculation of the bolometric light curve, 
and also showing the various individual radioactivity contributions. When we 
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Fig. 2. Bolometric light curve of SN 1987A with M( 56 Co) = 0.07 M & , 

M( 57 Co) = 3.3 x 10~ 3 Mq, and M( 44 Ti) = 1 x 10~ 4 M & . The separate contributions 
from the different isotopes as well as the positron input are shown as dashed and 
dotted lines. Note the dominant positron input from 44 Ti at late epochs. 

compare the bolometric luminosity with the total energy input one sees that 
between ~ 1200 days and ~ 2500 days the bolometric luminosity is up to a 
factor three higher than the input. At very late time most of the energy input 
is in the iron rich zones by positrons from 44 Ti, making time dependent effects 
less important. 

The most serious problem of using the bolometric light curve at late times is 
the large bolometric correction necessary because nearly all of the emission 
emerges in the unobserved far-IR. In figure 3 we nevertheless show the time- 
dependent bolometric light curve. This represents an extension of the modeling 
in Fransson & Kozrna (1993) to later epochs, and with updated physics. The 
freeze-out increases in this case the luminosity by a factor 2 - 3 at epochs 
later than 1200 days. The observations come from CTIO data by |5untzeff~et] 
al. (19921) , and ESO data by panziger fc Bouchet (19931) . While the agreement 
with the CTIO data is good, the agreement is considerably worse with the ESO 
data during the late epochs. The discrepancy between the ESO and CTIO data 
probably gives a good estimate of the systematic errors in the observations, 
and it is clear that a different approach has to be taken for estimating the 
mass of 44 Ti, and also to get a more accurate estimate of the 57 Ni mass. 

The second approach involves modeling of the line emission from the ejecta. 
An example of this was shown in figure 1, where we show Ha light curves for 
three different values of the 44 Ti mass. . Although the model reproduces the 
observations very well during the whole evolution, as a measure of the 44 Ti 
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Fig. 3. Bolometric light curve with M( 56 Co) = 0.07 M 0 , M( 57 Co) = 3.3 x 10 -3 Mq, 
and M( 44 Ti) = 1 x 10 -4 Mq. Data are from |Danziger fe Bouchet (1993| ) (dots) and 
from Suntzeff et al.(1992) (triangles). The dotted line gives the steady state light 
curve, while the solid shows the same with time dependent effects included. 


mass the light curve is clearly not very sensitive. Nearly all of the excitation 
is due to gamma-rays absorbed by the envelope, and the line is therefore 
strongly affected by the freeze-out. Most of the emission at < 3000 days is in 
fact due to the delayed recombination from the epoch at < 2000 days when 
57 Co dominated the energy input. A better alternative is to use lines from the 
metal rich core, in particular the Fe I and Fe 11 lines. An important result found 
by Chugai et ah (1997) is that in order to reproduce the Fe If flux, positrons 
have to be trapped efficiently in the iron rich gas. Pure Coulomb collisions 
are not sufficient for this, but a weak magnetic held is enough to make the 
gyro radius of the positrons small enough. A severe problem is, however, that 
modeling of individual lines is very sensitive to atomic data, which in spite of 
dramatic improvements from the IRON-project, is still a problem. In a similar 
attempt [Lundqvist et ah (2001| ) have tried to use the FIR [Fe II] 26.0/i fine 
structure line in connection with observations by ISO. Unfortunately, only an 
upper limit to the flux was obtained, yielding an upper limit to the 44 Ti mass 
of 1.1 x 10 ” 4 M 0 . There is, however, a number of loop holes in this result, such 
as extra dust cooling, mixing with other heavy elements contributing to the 
cooling like silicon and sulphur, as well as uncertainties in the atomic data. 


A compromise between these two approaches is offered by broad band photom¬ 
etry, where photometry in the optical and near-IR is compared to synthetic 
photometry from spectral calculations. This method avoids the need to ex¬ 
trapolate into the far-IR for the bolometric light curve, as well as many of the 
uncertainties involved in the calculation of individual lines. Rather than mod¬ 
eling individual lines, the broad band fluxes are mainly sensitive to the division 
between UV, optical and IR in the spectral distributions of the elements. This 
is determined by the atomic configuration, i.e. the term diagram, and less 
by the exact path a recombining photon takes when cascading to the ground 
state. Except for the fraction of the emission in the optical and IR which is 
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Fig. 4. B and V light curves from Suntzeff et al.(1991, 2001), together with models 
with (0.5,1.0, 2.0) x 10” 4 M 0 of 44 Ti. All models have 3.3 x 10~ 3 M 0 of 57 Ni 

absorbed by the dust, the photometry is insensitive to this component. The 
fraction absorbed by dust can be estimated from mid-IR observations, as well 
as from line profile shifts (Lucy et al. 1991). 

In figure 4 we show the B and V band light curves together with observa¬ 
tions. The reason for concentrating on these bands is that they are strongly 
dominated by lines from Fe I and Fe II, and therefore most sensitive to 
the 44 Ti mass, in contrast to the R band which is dominated by Ha, and 
therefore sensitive to the freeze out. Observations are taken from Suntzeff 
et al.(1991) and HST observations from Suntzeff et al.(2001). The models 
have 6.9 x 10 -2 M 0 of 56 Co, 3.3 x 10 -3 M 0 of 5 'Co and three different val¬ 
ues of the 44 Ti mass, (0.5,1.0, 2.0) x 1CT 4 M 0 . As can be seen, both the B 
and V bands are well reproduced by the observations throughout the whole 
evolution, and the currently best estimate of the 44 Ti mass is in the range 
(0.5 — 2.0) x 10 -4 Mq. This method also yields the best estimate of the 
57 Ni mass, ~ 3.3 x 10 -3 M & . A more complete discussion of these results 
will be given in Fransson & Kozrna (2001). 


4 Conclusions 

SN 1987A is by far the best source of nucleosynthesis of any SN. Reliable 
masses for both the most abundant elements, H, He, and O, as well as the the 
most abundant radioactive isotopes, 56 Ni, 57 Ni, and 44 Ti have been obtained. 
These provide some of the most important constraints on the progenitor and 
on the explosion itself. 

CF is grateful to Roland Diehl for organizing a very successful workshop, and 
to Don Clayton and Nikos Prantzos for useful discussions and comments. 
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